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PREFACE 


This  Memorandum  is  one  result  of  a  continuing  investigation  of 
the  radar  returns  from  missile  wakes.  In  a  companion  Memorandum, 
RM-3573-ARPA,  The  Scattering  of  Electromagnetic  Waves  from  Plasma 
Cylinders:  Part  I.  a  general  theory  was  developed  for  the  short- 
wavelength  limit  of  the  electromagnetic  radiation  by  radially  vary¬ 
ing  electron-density  distributions.  This  Memorandum  considers  the 
techniques  for  calculating  the  long -wave length  limit  of  the  scatter¬ 
ing. 

The  work  is  part  of  RAND's  study  of  midcourse  phenomenology, 
conducted  under  the  Advanced  Research  Projects  Agency's  Defender 


program. 
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SUMMARY 


In  this  Memorandum  we  consider  the  long -wavelength  limit  of  the 
scattering  of  electromagnetic  radiation  in  the  transverse  magnetic 
mode  by  infinite  plasma  cylinders  with  radially  varying  electron- 
density  distributions-  Techniques  for  calculating  the  scattering 
cross  section  are  described  for  both  increasing  and  decreasing 
electron-density  distributions.  Resonance  phenomena  peculiar  to 
increasing  electron-density  distributions  are  discussed  in  some 
detail  with  specific  numerical  examples. 

Finally,  we  present  a  resume  of  methods  one  can  use  to  calcu¬ 
late  the  scattering  cross  section  for  radiation  of  arbitrary  wave¬ 
length  for  monotonically  decreasing  elecltr on-density  distributions. 
Examples  are  given  for  the  Gaussian  electron  distribution. 
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I.  INTRODUCTION 

•ff 

In  a  previous  Memorandum,  hereafter  referred  to  as  Part  I,  we 
developed  a  general  theory  for  the  scattering  of  electromagnetic  radia¬ 
tion  by  radially  varying  electron-density  distributions  under  the 
conditions  that  the  incident  wave  was  a  plane  wave  traveling  perpen¬ 
dicular  to  the  cylinder  axis  and  polarized  parallel  to  the  cylinder 
axis,  the  electron-density  distribution  was  monotonic,  and  collisions 
were  negligible.  Techniques  for  calculating  the  cross  section  in 
the  short-wavelength  limit  were  described  and  applied  to  a  few 
specific  cases. 

In  this  Memorandum  we  will  consider  the  long -wave length  limit 
of  the  scattering.  The  conditions  of  the  scattering  and,  therefore, 
the  formulation  of  the  problem  are  the  same  as  in  Part  I.  However, 
the  techniques  appropriate  for  calculating  the  scattering  cross 
section  differ  in  the  long-  and  short-wavelength  limits- 

In  our  consideration  of  the  long-wavelength  limit,  we  will 
place  special  emphasis  on  the  scattering  by  monotonically  increasing 
electron-density  distributions  (analogous  to  the  scattering  by 
attractive  potentials  in  particle-scattering  theory),  since  distri¬ 
butions  of  this  type  have  not  been  previously  considered  in  the 
literature  except  in  the  small  perturbation  limit.  Such  distribu¬ 
tions  can  give  rise  to  anomalies  in  the  scattering  not  encountered 
with  decreasing  electron-density  distributions.  These  anomalies 
are  associated  with  partial-wave  resonances,  familiar  in  the 

* 

Phyllis  Grei finger,  The  Scattering  of  Electromagnetic  Waves 
From  Plasma  Cylinders:  Part  I.  The  RAND  Corporation,  RM-3573-ARPA, 
August  1963. 
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analogous  situation  of  the  scattering  of  low-energy  particles  from 
spherically  symmetric  attractive  potentials. 

Despite  the  emphasis  placed  on  the  increasing  electron- density 
distributions,  the  methods  described  for  calculating  the  cross 
section  in  the  long -wave length  limit  will  be  equally  applicable  to 
increasing  and  decreasing  electron-density  distributions. 
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II.  REVIEW  OF  PARTIAL-WAVE  EXPANSION 


The  formulation  of  the  problem,  the  expansion  of  the  wave  into 
partial  waves,  and  the  expression  of  the  scattering  cross  section  in 
terms  of  the  phase  shifts  of  the  partial  waves  were  presented  in  some 
detail  in  Part  I.  We  will  present  a  brief  review  here,  with  a  few 
changes  in  notation  which  are  introduced  for  convenience  in  consider¬ 
ing  the  long -wave length  limit  of  the  scattering. 

_itut 

Dropping  the  time-dependent  factor  e  ,  the  wave  equation  in 
cylindrical  coordinates  P,  Z,  and  cp  for  the  electric  field  E  ■  E^  is 


1  1SE 

,2  P  c>p 


(u,2  - 


plasma  frequency 


Letting 

y  -  P/P0 
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Eq.  (1)  becomes 


52E  1  3E  l__  a^E 

by2  y  Zy  y2  a?2 


E  -  0 


(2) 


The  characteristic  radius  of  the  scattering  region  is  Pq,  and  the 
subscripts  0  and  ®  on  ng  and  U)p  refer  to  p  ■  0  and  p  -  °°,  respectively. 
The  dimensionless  parameters  K  and  3  are  related  to  the  wave  number 
k  and  the  parameter  or  of  Part  I  by 


K  -  kp 

o 

(3) 

2  2  2 
3  «  k  p  Of  -  K  Of 
o 


The  parameter  3  is  not  a  function  of  the  frequency  of  the  radiation. 

For  monotonic  electron- density  distributions,  it  depends  upon  the 
difference  between  the  maximum  and  minimum  electron  densities  as  well 
as  on  the  size  of  the  scattering  region.  It  is  a  more  convenient 
parameter  than  for  small  K,  since  it  remains  finite  as  K  "*  0.  The 
small-K  limit  refers  to  /u)j2(u>Po/c  j  <  1.  Note  that  this 

situation  can  occur  even  if  «jpQ/c  is  large,  provided  that  <d  Is  not 
much  larger  than  u>  •  For  example,  for  the  ionosphere,  is  of  the 

-f-Q  J  2 

order  of  10  rad/sec  (ng  ~  3  x  10  /cm  ) .  Large  negative  values 
of  3  can  occur  for  electron- density  "holes"  such  as  the  partially 
evacuated  region  behind  a  hypersonic  body  moving  through  the  ionosphere 
if  the  body  radius  is  larger  than  a  few  meters.  For  K  to  be  less  than 
one  under  these  conditions,  we  require  radiation  very  close  in  fre¬ 
quency  to  the  ionospheric  plasma  frequency,  l.e.,  of  the  order  of 
20  Mc/sec. 


-5- 


Expanding  E  in  a  Fourier  series,  we  get 


E  -  E  c  cos  (tncp)  *  (y) 
m-0  m  “ 

C  *  1  m  ■  0  (4) 

m 

e  ■  2  m  ^  0 
m 


The  function  ♦  (y),  which  is  related  to  u  (y)  of  Part  I  by 
,  »  m  m 

um(y) 

»m(y)  *  ~^aw=~,  satisfies  the  differential  equation 


♦"(y)  +  -  ♦  *Cy)  + 

“  J  y  m w 


m 


K2  -  Pg(y)  -  ^ 

y  J 


♦  (y) 

m 


(5) 


with  the  boundary  conditions 


♦  (o)  is  finite 
m 


ms 


♦  -  If  .  constant  x 

m  mi 


elK, 

7T 


as  y 


The  subscripts  8  and  i  refer  to  the  scattered  and  incident  wave, 
respectively.  If  the  asymptotic  form  of  t  is  written  as 

i  _/  2  \1'2  imTT/2  -IH  /„  mn  n  -  \  ,,, 

(ml)  e  e  “  CO*  (Ky  -  —  -  4  -  \j  (6) 

where  T]  is  the  phase  shift  of  the  m^  partial  wave,  the  scattered 
m 

field  E(  for  large  y  becomes 
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ei(Ky-TT/4)  ®  /  -21“^ 

'  2  «  cos  raSple 

yZTiKy  m«0  m  T  \ 
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eiKy 
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(7) 


The  differential-scattering  cross 


l2 

section  f (cp)  J  then  is  given  by 


»  »  o  P 

f=" 

T-l 

CM 

f 

8 

lf(cH  “  2^K 

m-0  *m  COS  (nKp)\e  “  “  VI 

(8) 


It  was  shown  in  Section  IV  of  Part  I  that  the  smaller  the  K, 
the  fewer  the  phase  shifts  which  contribute  to  the  scattering.  As 
K  "♦  0,  we  can  generally  expect  all  the  phase  shifts  except  Tj  to  be 
negligible.  The  resultant  scattering  will  be  isotropic,  with  |f(<p)|^ 
given  by 


|fOrt|2^»ln2  1). 


(9) 


As  we  will  see,  exceptions  can  arise  when  it  is  necessary  to  consider 
one  or  more  higher  phase  shifts,  even  as  K  0.  The  exceptions  occur 
when  either  |  T)o  |  is  an  integral  multiple  of  n,  causing  sin  Tj^  to 
vanish,  or  when  one  of  the  higher  m  partial  waves  is  in  resonance, 
causing  a  phase  shift  which  is  an  odd  multiple  of  tt/2. 

If  K  is  small  but  finite  and  we  do  not  have  a  resonance  situation 
for  m  >  1,  we  need  consider  only  and  1)^.  The  differential  cross 
section  I  f  Cp)  1 2  becomes 


|f(cp)j  ^  —  fjjjr-  Isin^  1)o  +  4  cos^  cp  sin^  +  4  cos  cp  (sin^  Tlo  sin^  T|^ 

\ co*  \)] 


(10) 


+  sin  11  cos  T)  sin 
'o  o 


-7- 


We  can  see  from  Eq.  (10)  that  the  ratio  of  the  contribution  of  the 
m  -  1  term  to  the  m  -  0  term  to  the  total  cross  section 

2TT 

J*  I  f  (cp)  I  2  dtp 

o 

2 

is  proportional  to  sin  T]^.  However,  the  ratio  of  the  largest  angle- 
dependent  term  to  the  constant  term  in  the  differential  cross  section 
is  proportional  to  sin  Thus,  the  partial  wave  m  -  1  generally 
manifests  itself  in  the  differential  cross  section  at  a  smaller  value 
of  K  than  that  at  which  it  becomes  significant  in  the  total  cross 
section.  For  example,  if  T|  »  45  deg  and  T|^  -  4-5  deg,  the  m  ■  1 
wave  contributes  only  about  2.4  per  cent  to  the  total  scattering, 
while  it  makes  the  forward  scattering  (<p  ■  0)  about  twice  as  large 
as  the  backward  scattering  (cp  «  tt)  . 
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III.  SMALL-K  PHASE  SHIFTS 

FINITE-RADIUS  CYLINDER  SOLUTIONS 

Let  us  consider  functions  g(y)  which  vanish  for  y  >  1.  For 
such  cases,  the  Eq.  (5)  solutions  inside  and  outside  the  cylinder 
are  matched  at  the  boundary;  for  the  T.M.  mode,  both  (y)  and 
^(y)  are  continuous  at  y  =  1.  The  phase  shifts  are  then  given 
by 


cot 


KH’(K)  -  VH(I) 
tu  mm 

V  J  (K)  -  KJ'OO 
m  m  7  m 


m 


(ID 


The  functions  J^CK)  and  N^K)  are  Bessel  and  Neumann  functions,  respec¬ 
tively.  The  Ym's  are  found  from  the  solutions  to  Eq.  (5)  inside  the  cylin¬ 
der;  for  each  m,  the  required  solution  is  that  which  is  finite  at  y  -  0. 

For  K  «  1,  we  have  approximately 


cot  T|  a:  1 
o  tt 


C  -  Euler ' s  constant 


-(,  +  YjHm(K)  +  KViW 


cot  1^  = 


m*.  (m  -  1) 
n 


(12a) 


(12b) 

(12c) 


(12d) 


-9- 


We  will  discuss  two  cases  £or  which  we  can  obtain  analytic  expressions 

for  the  V  1  s. 
m 

Homogeneous  Cylinder 

Equation  (5)  can  be  solved  analytically  for  a  step-function 
electron-density  distribution  (homogeneous  cylinder) : 


g(y)  “1  y  *  1 

=  0  y  >  1 

The  solutions  are  the  Bessel  functions  J  (Jv?  -  P  y)  regular  at  the 
origin.  The  Vm'8  are  given  by 


m 


(13) 


Using  the  recurrence  relationships  of  the  Bessel  functions,  we 
get  for  m  2  1 


Y  +  m 
m 

Y  -  m 
in 


j  XU2  -  p) 


mfl' 


(Vk2  -  p) 


(14) 


so  that  an  approximate  expression  for  cot  for  K  «  1  and  m  *  1  Is 


2m  J 


-  ») 


(15) 


2  2 
For  a  given  small  K  and  P  <  K  ,  there  are  maxima  and  minima  for  sin 

(m  s  1)  which  occur  at  the  zeroes  of  J  .(Vk2  -  3 )  and  J  ,.(Yk2  -  p), 

m-  i  m+l  * 


respectively. 
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Quadratlc  Distribution 

Another  electron-density  distribution  for  which  we  can  obtain  an 
analytic  solution  to  Eq.  (5)  is  the  so-called  quadratic  distribution 

g(y)  -  1  -  y2  y  S1 

-  0  y  ^  1 


The  solution  inside  the  cylinder  for  this  case  is 


i Ji 


■  e 


2 

L- 

2 


*Vi(v  v  •  ‘J* y2) 


(16) 


m  +  1  i  U2  -  p) 
a»-  — 


c  »*  +  1 
m 


The  function  jF^(a;  b;  Z)  is  the  confluent  hypergeometrlc  function 


^Ui  b;  Z) 


1 


a  Z 

+  bIT  + 


a(a  +  1)  Z2  , 

b(b  +  1)  IT  + 


(17) 


The  Ym's  can  be  found  from  the  recursion  formulas  for  the  hyper- 
geometric  functions  and  are  given  by 


-  We 


i*i  (y  H 

i*i(m«:  v  -  &)  . 


(18) 


EVALUATION  OF  T)  FOR  ARBITRARY  ELECTRON'* DENSITY  DISTRIBUTIONS 
o 

With  certain  exceptions  to  be  discussed  later ,  Eq.  (9)  repre¬ 
sents  e  fair  approximation  to  the  small-K  cross  section.  It  would 
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be  especially  useful,  therefore,  to  develop  a  technique  for  calcu¬ 
lating  T\q  for  arbitrary  functions  g(y)  (subject  to  the  restriction 

2 

that  g(y)  falls  off  faster  than  1/y  for  large  y) •  We  will  do  this 
by  analogy  with  a  technique  used  for  obtaining  an  approximation  to 
the  s-wave  phase  shift  for  the  scattering  from  spherically  symmetric 
scattering  regions* 

Let  us  rewrite  Eq.  (5)  for  m  ■  0  as  follows 

&(y-af)+(K2  -  Bg(j,))y,oK  ■ 0  <w> 

with  boundary  conditions 


♦^(0)  finite 


Vy> 


cos  (ity  -  J  - 


as  y  -*  <*> 


The  subscript  K  simply  refers  to  the  solution  for  a  particular  K. 

The  function  toQ(y)  would  be  the  solution  to  Eq.  (19)  (finite  at  the 
origin)  with  K  -  0. 

Let  us  also  consider  functions  W^(y)  which  have  the  same 
asymptotic  behavior  as  toR(y)  and  which  satisfy  Eq.  (19)  with  3*0 

%(y^)+K2yWK‘° 


Vy)  -  Vy) 


as  y  -*  • 


(20) 
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Equation  (28)  can  be  solved  immediately,  and  we  get 

W  (y)  -  B  fcot  T)  J  (Ky)  +  N  (Ky)l  for  K  4  0 
K  L  o  o  o  J 


ty 

The  factor  —  is  put  in  the  function  W  (y)  to  cause  yW’(y)  and  yW'(y) 

TT  O  O  r^. 

to  be  equal  at  y  ■  0.  The  parameter  V  can  only  be  found  from  the 

1  *o 

solution  of  Eq.  (19)  with  K.  -  0.  Since  V  is  equal  to  ~  -jr—  at  y  -  1, 

o  y 

it  is  the  same  as  VQ  of  Eq.  (11)  for  K  =  0  and  for  those  functions 
g(y)  which  vanish  for  y  >  1.  The  amplitude  B  is  as  yet  unspecified, 
since  we  have  not  specified  the  value  of  at  y  =  0.  We  will 
choose  B  by  requiring,  for  convenience,  that  i|fQK(0)  *  1.  From  the 
differential  Eq.  (19),  we  then  find  that  [ ^(y) /y ]  ■  P/2-  Obtaining 
numerical  solutions  of  Eq.  (19)  with  K  -  0  is  a  straightforward 
problem  for  functions  g(y)  which  fall  off  faster  than  for  large  y. 

y 

From  the  solutions,  we  obtain  the  parameters  B  and  V  as  functions  of 
0- 

Let  us  now  see  how  we  can  obtain  cot  T)  for  arbitrary  small  K 
from  the  function  \J/qoC y)  -  Let  us  consider  the  following  two  equalities 


Subtracting  Eq.  (22b)  from  Eq.  (22a),  integrating  over  y  from  0  to  ®, 
and  inserting  the  boundary  conditions,  we  get 


(23) 


=  £*2 


2 


■  r2  r(«A  -  ** 


From  this  exact  equation,  an  approximate  expression  for  cot  “H 

2 

correct  to  first  order  in  K  may  be  obtained  by  replacing  Wv  by  W 

K  i 

and  „  by  t  in  the  integral 
oK  J  oo  ° 


cot  T) 

o 


K2A 


-?r 


y  +  i)2  - 


♦oo  <T> 


B 


y  dy 


(24) 


Equation  (24)  can  be  expected  to  break  down  as 

Y  —  0  (|  cot  Tl^  J  -*  oc^ 

For  values  of  0  which  are  such  that  we  are  near  a  zero  of  y,  we 

cannot  replace  the  wave  functions  for  finite  K  by  those  for  K  ■  0, 

since  small  changes  in  K  have  large  effects  on  the  wave  functions. 

The  zeros  of  sin  Tj  for  small  K  will  cause  the  cross  section 
o 

to  exhibit  deep  minima.  If  0  is  negative,  the  smallest  value  of 

|0|  for  which  sin  T]  ■  0  will  be  that  which  will  cause  the  m  ■  0 
'  ■  o 

wave  to  be  pulled  in  by  just  half  a  cycle,  giving  a  phase  shift 
T|  of  -tt  as  illustrated  in  the  following  sketch. 
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The  effect  is  a  diffraction  of  the  wave  around  the  cylinder, 
in  which  the  wave  inside  the  cylinder  is  distorted  in  such  a  way 
that  it  fits  on  smoothly  to  an  undistorted  wave  outside.  This  par¬ 
ticular  phenomenon  with  K  small  is  associated  only  with  negative  0; 
for  positive  0,  T]q  is  n  (m  »  0  wave  pulled  out  by  half  a  cycle) 
only  for  fairly  large  K,  so  that  the  higher  m  phase  shifts  are  not 
negligible.  Such  minima  in  the  cross  section  are  well  known  in  the 
analogous  situation  of  the  scattering  of  low-energy  particles  from 
spherically  symmetric  attractive  potentials.  In  particular,  an 
extremely  low  minimum  known  as  the  Ramsauer -Townsend  effect  occurs 
in  the  scattering  cross  section  of  low-energy  electrons  by  rare-gas 
atoms. 

When  we  are  not  at  such  a  minimum,  we  can  see  from  Eq.  (24) 
that  for  finite  P,  as  K  -*  0 

cot  l|o  -  |  tv(±)  (25) 


so  that 
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?  np„  /  l 

lf«P)|  -  ar(^  £)  (26) 

For  small  but  finite  K,  it  is  necessary  to  use  the  complete  Eq.  (24) 

to  evaluate  cot  1]  . 

o 

The  parameter  Y  and  the  integral  A  were  computed  on  the  IBM  7090 

_  2 

for  a  Gaussian  distribution  g(y)  *  e  ^  for  a  variety  of  values  of 

1  2 

positive  and  negative  0.  In  Fig.  1,  we  have  shown  —  and  —  A  versus  0 
for  the  Gaussian  distribution  for  positive  0.  In  Fig.  2,  we  have 
plotted  ^  versus  0  for  negative  0  for  the  Gaussian  distribution, 
homogeneous  cylinder,  and  quadratic  distribution  for  those  values 
of  0  for  which  |^|  <  2.  In  Fig.  3,  we  have  plotted  J^j  versus  0  for 
the  homogeneous  cylinder  and  Gaussian  distribucion  in  the  region  in 
which  Y  -*  0.  The  integral  A  for  the  Gaussian  distribution  is  shown 
as  a  function  of  0  for  negative  0  in  Fig.  4. 

The  approximate  cot  TjQ,  given  by  Eq.  (24),  and  the  exact  cot  T)o 
were  calculated  for  a  homogeneous  cylinder  for  0  *  -4,  -14.2,  and 
-14-6.  For  0  *  "^(y  “  -0.194),  it  was  found  that  the  error  involved 
in  the  approximation  was  less  than  1  per  cent  for  K  up  to  unity,  and 
the  results  were  not  plotted.  For  0  ■  -14.2  and  -14.6,  |y|  is  small, 
and  one  can  expect  the  approximation  given  by  Eq.  (24)  to  break  down 
as  K  -»  1.  In  Fig.  5,  we  have  plotted  both  the  approximate  and  exact 
values  of  cot  T|q  versus  K  for  0  ■  -14-2  and  -14.6.  For  0  •  -14.2^  ~  5) , 
the  approximation  is  surprisingly  good  for  K  near  unity;  the  error 
in  cot  T)^  is  only  20  per  cent  for  K  ■  1.  For  0  -  -14.6^  «  22.3), 
the  approximation  breaks  down  badly  as  K  -*  1  but  is,  nevertheless, 
excellent  for  small  K;  for  K  *  0.5,  the  error  is  only  8.5  per  cent. 
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Fig.  1— Curves  of  —  and  A  versus  P  for  p  positive 


20 


-17- 


Fig.  2— Curves  of  y  versus  (3  for  negative  (3  ( |  y  |  <  2) 


-18- 


-P 


Fig.  3— Curves  of  |y|  versus  p  for  negative  p  near  y  0 
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Fig.  4 — Curves  of  Jr  A  versus  £  for  negative  £ 
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Fig.  5— Comparison  of  exact  and  approximate  values 
of  cot  rj0  for  lyl  large 
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IV.  RESONANCE  SCATTERING 

We  can  see  from  Eq.  (11)  for  the  finite-radius  cylinder  that 
N^(K) 

when  Y  “  K  ~  .  cot  T|  will  be  zero.  Whenever  cot  Tj  ■  0,  the 

m  N  (K.)  m  m 

th  m 

m  partial  wave  is  said  to  be  in  resonance. 

Small-K  resonances  for  m  s  1  cannot  arise  if  @  is  positive.  We 

can  easily  prove  this  for  finite-radius  cylinders  by  showing  that 

N^(K) 

for  @  >  0,  m  ^  1,  and  K  <  1,  y  is  always  positive.  Since  ^  -r-r  <  0 

m  m'  ' 

for  small  K,  a  resonance  can  only  occur  if  Ym  is  negative. 

Let  us  rewrite  Eq.  (5)  as  follows 
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infinite  sequence  of  values  0  of  negative  0  for  which  Y  will  have 

mn  in 

precisely  the  right  negative  value  to  cause  T|  ■  -(2n  -  1)  where 

n  is  a  positive  integer.  For  example,  for  the  homogeneous  cylinder, 

the  0  's  are  found  from  the  roots  of  the  equation  J  -(/ K2  -  0  )  ■  0. 

mn  m-lx*  on/ 

n - 

Considered  as  a  function  of  m,  the  lowest  root^K  -  0^^  will  increase 
with  m. 

In  a  small-K  resonance  situation  .for  m  ^  1,  the  scattering  will 
usually  be  the  result  of  the  interference  between  the  wave  in  resonance 
and  the  m  *  0  wave.  The  differential  cross  section  in  the  vicinity 
of  an  i  ^  1  resonance  will  be  approximately 

2p  I" 

|f(cp)|2  —  |sin2  T)q  +  4  cos2  «p  sin2 


+  4  cos  n*p(sin2  T)o  sin  +  sin  sin  cos  T)^  cos 


v] 


(29) 


The  angular  distribution  S(cp)  at  the  peak  of  the  resonance  will  be 
given  by 


S(cp) 


\  Usl  \ 


Sf.TVwi  2 


2  2  2 
sin  T1  +  4  cos  mcp  +  4  cos  oq>  sin  T)o 

sin2  T1  +  2 
o 


(30) 


Illustrations  of  resonant  scattering  are  shown  in  Figs.  6  to  9. 

In  Fig.  6,  we  have  shown  |  sin  Tj^j  and  |sin  T)o|  as  functions  of  K  near 

an  m  ■  1  resonance  for  a  homogeneous  cylinder.  In  Fig.  7,  we  have 

plotted  the  radar  cross  section  |  f (tt)  |  near  the  resonance  in  com- 

2P  2 

par  is  on  with  the  small-K  cross  section  sin  f)o  in  the  absence  of 
ana  2  1  resonance.  Figure  8  shows  the  angular  distribution  S(cp)  at 

2 

the  peak  of  this  particular  m  ■  1  resonance  and  compares  it  to  2  cos  cp. 
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the  angular  distribution  of  a  pure  m  ■  1  wave.  In  Fig.  9,  we  have 

shown  the  angular  distribution  at  the  peak  of  an  m  -  1  resonance  for 

2 

the  quadratic  distribution  in  comparison  with  the  2  cos  <p  distri¬ 
bution.  Figures  8  and  9  show  the  effects  of  the  interference  between 

the  m  »  0  and  m  -  1  waves  on  the  angular  distribution.  The  distor- 
2 

tion  of  the  2  cos  cp  distribution  by  the  m  =  0  wave  is  more  pro¬ 
nounced  in  Fig.  9  than  in  Fig.  8  simply  because  |  sin  T)  |  is  larger 
in  the  former  case. 

In  Fig.  10,  we  have  shown  |  sin  T|q  |,  j  sin  T]^  j,  and  j  sin  T|^  | 
near  an  m  -  2  resonance  for  the  homogeneous  cylinder.  It  happens 
that  this  particular  resonance  occurs  near  a  minimum  of  |  sin  Ho|, 
so  that  at  the  peak  of  this  m  *  2  resonance  we  have  an  almost  pure 
m  -  2  wave.  This  result  is  not  surprising.  We  can  see  from  Eq.  (15) 
that  for  the  homogeneous  cylinder,  cot  Tj^  sa  0  when  -i/Vk2  -"  b)«  0. 
However,  from  Eq.  (12)  we  see  that  for  small  K,  |  cot  T|o|  -*  ®  when 

J*2  -  e  jX/k2  -  @) 

is  very  small  in  magnitude  or  when  J  i-Jxf  -  is  near  zero.  We  have 

not  plotted  the  angular  distribution  at  the  peak  of  the  particular 

m  ■  2  resonance  being  considered  here,  since  it  is  essentially  a  pure 
2 

2  cos  29  distribution. 

So  far  we  have  discussed  only  resonance  maxima  for  the  m  S:  1 

partial  waves.  We  may  also  have  smallrK  resonances  for  m  ■  0;  they 

1  2 

occur  in  the  vicinity  of  r  ■  -  in  — r  (see  Eq.  (24))..  If  K  is  very 

Y  KeC 

small,  the  m  -  0  resonance  maxima  will  occur  for  y  negative  and  fairly 
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2 

small  in  magnitude;  however,  we  have  seen  that  the  minima  in  sin  T| 

occur  for  |y|  very  small*  Therefore,  for  a  given  small  K,  if  we 
2 

consider  sin  T]  as  a  function  of  &  for  negative  (3  in  a  region  where 

2 

Y  passes  through  zero,  we  can  expect  a  zero  and  a  maximum  in  sin  Mo 

to  be  fairly  close  together.  We  have  shown  this  effect  in  Fig.  11  by 
2 

plotting  sin  H0  versus  0  for  a  homogeneous  cylinder  for  K  »  0.06  and 
-0  between  13  and  16.  We  can  see  that  a  minimum  occurs  for  -0  near 
14.7  and  a  maximum  for  -0  near  15.4.  For  K  smaller,  the  minimum  and 
maximum  would  be  still  closer  together. 
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Fig.  6— Curves  of  m  •  Oand  m  ■  1  phase  shifts  near  an  m  *  1  resonance 
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Fig.  7—  Effect  of  an  m  -  1  resonance  on  the  radar  cross  section 
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4> 


Fig.  8— Angular  distribution  at  the  peak  of  an  m  •  1  resonance 
(homogeneous  cylinder;  /?  -  -5.76,  K  -  0.058) 
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Fig.  9— Angular  distribution  at  the  peak  of  an  m  -  1  resonance 
(quadratic  distribution;  p  ■  -12,  K  -  0.25) 
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Fig.  10— Phase  shifts  near  an  m  -  2  resonance 
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V.  CROSS  SECTION  FOR  I  3 1  <  1  AND  K  <  1 

Since  |3|  and  K  are  both  proportional  to  pQ,  in  many  situations 
in  which  K  is  small,  |3l  will  also  be  less  than  one.  It  would  be 
useful,  therefore,  to  develop  a  simple  method  for  calculating  the 
cross  section  for  |3|  <1  and  K  <  1. 

We  know  that  for  K  finite,  we  can  use  the  Born  approximation  to 
estimate  the  cross  section  or  the  individual  phase  shifts  for  |3l  -*  0. 
For  example,  the  Born  approximation  for  "H  is  given  by 

yg(y)Jo(Ry)  dy  (3i) 

o 


If  we  let  K  go  to  zero,  this  gives 


\  i  e  *1*°  y s(y)  dy 


(32) 


However,  we  have  seen  that  for  3  finite,  if  we  let  K 
so  that 


°»  cot  \  n  ln  b 


tt 


°  2ln* 


(33) 


Therefore,  for  K  -*  0,  we  cannot  use  Eq.  (32)  for  finite  3,  even  for 

I  3  I  small.  We  need  a  simple  approximation  to  Tj^  for  K  <  1  and 

|3l  <  1,  which  reduces  to  Eq.  (32)  for  K  finite  and  |3|  -»  0,  and  to 

Eq.  (33)  for  I  3  I  finite  and  K  -*  0. 

Since  Eq.  (24)  represents  a  good  approximation  to  T|  for  small 

o 

K  and  all  finite  values  of  |  3 j  (except  at  the  cross-sectional  minima), 
we  can  obtain  the  approximation  we  are  seeking  by  obtaining  a  power 
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series  expansion  in  0  of  1/Y  and  A.  This  can  be  done  by  expanding 

to  (y)  in  a  power  series  in  0 
oo 


*  nio  6\<y) 


(34) 


From  the  differential  equation  for  t00>  we  get  the  equations  for 

the  functions  Y  (y) 
n 

%(y£)-o 

%(y  df)  ■  y*<y>Y„-l<y>  " 2 1 

Y  (0)  *  1 
o 

Y  (0)  -  0  n  0 

n' 

Y'  (0)  *  0 
n 


The  solutions  then  are 

Vy>  ■ 1 


y  (y)  ■  7  y's(y’)  4n  y  1(y')  dy'  (36) 

o  y 

n  s  1 

We  can  see  from  Eq.  (36)  that  for  n  ^  1 

Y  (y)  -  C  In  y  +  D 
nv  1  n  J  n 


as  y  -*  * 


(37) 
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where 


cn  =  f°  y,g(y,)Yn_1(y')  dy' 


Dn  =  -  J°°  y'g(y')  y'  Yn  l(y')  dy' 


Since  t  (y) 
oo w 


2  1 

—  B  (4n  y  +  — )  as  y  -*  00 ,  we  have 


|  B  -  E  0nC 
n  n*l  n 


(38) 


|  ~  =  i  +  e  e”D 

tt  V  n*l  n 


1 

Y 


1  +  £,8*0 

n*l  n 

oo 

E.  8nC 
n=l  n 


We  see  that  P/Y  -*  —  as  8  -♦  0.  The  next  approximation  yields 
C1 


(39) 


The  integral  A  in  Eq.  (24)  is  proportional  to  g-  for  small  B,  and  the 
leading  term  is 


a  -  -  — o  r  y3g(y)  dy 

tt8c  j  o 


(40) 
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Substituting  Eqs.  (39)  and  (40)  into  Eq.  (24),  we  get 


cot 


+  -^~2  r  y3s(y)  dy 

20C*  o 


(41) 


If  we  let  0  remain  fixed,  then  as  K  -•  0,  Eq.  (41)  reduces  to  Eq.  (33). 
For  K  small  but  finite,  as  P  ■*  0,  Eq.  (41)  becomes 

\  "*  f  ^l^1  ■fc^J'00  y3g(y)  dy)  "  f  P[j"  yg(y)  dy  “  r*  y3gCy)  dyj 

2 

This  is  the  s->rae  as  Eq.  (31)  if  we  expand  JQ(Ky)  and  keep  the  first 
two  terms. 

Thus,  we  can  expect  Eq.  (41)  to  represent  a  good  approximation 
to  T)  for  K  <  1  and  <  1  for  those  functions  g(y)  which  fall  off 
rapidly  enough  at  large  distances  to  insure  convergence  of  the  inte¬ 
grals. 

Let  us  note  from  Eq.  (41)  that  for  negative  0  we  can  have  an 
m  *  0  resonance  even  for  1$)  small.  For  a  resonance  to  occur  for 
|p|  <  1,  K  is  given  by 


(42) 


In  Fig.  12,  we  have  shown  the  correct  radar  cross  section  for 
K  ■  0.06  for  a  homogeneous  cylinder  for  small  0  and  compared  it  to 
the  Born  approximation.  The  peak  at  -0  ~  0.63  is  due  to  an  m  ■  0 
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resonance,  and  it  can  be  seen  that  the  Born  approximation  is  very 
poor  near  this  resonance.  The  fact  that  for  the  particular  case 
considered  the  Born  approximation  cross  section  happens  to  be  good 
at  3  ■  -1  is  fortuitous. 
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Fiq.  12 — A  comparison  of  the  Born  approximation  with  the  correct 
radar  cross  section  for  small  1/31  and  K 
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VI .  Rffsiftlf 

In  Part  I,  we  described  the  Born  approximation  and  the  geometrical 
optics  approximation  to  the  scattering  cross  section  and  discussed 
their  limits  of  applicability.  We  also  described  methods  for  approxi¬ 
mating  the  individual  phase  shifts  for  those  values  of  large  K  for 
which  the  Born  or  geometrical-optics  approximations  are  not  applicable. 
In  this  Memorandum  we  have  confined  ourselves  to  a  consideration  of 
the  cross  section  for  K  <  1. 

The  following  briefly  summarizes  the  methods  one  can  use  to 
calculate  the  cross  section  for  different  K  values  for  |0|  <  1  and 
0  >  1.  Since  we  did  not  consider  specific  methods  for  calculating 
the  cross  section  for  0  <  -1  and  K  >  1  in  either  Part  I  or  Part  II 
(except  in  the  Born -approximation  limit),  we  will  not  include  0  <  -1 
in  this  summary. 

lei  <  i 

For  small  |0|  and  K  <  1,  we  may  use  Eq.  (24)  for  the  evaluation 
of  cot  1]  and  the  Born  approximation  to  the  higher  phase  shifts- 
For  larger  values  of  K,  the  Born  approximation  to  the  differential 
cross  section  is  applicable.  Thus,  for  |0|  <1,  we  can  get  a  good 
estimate  of  the  cross  section  for  all  K. 

As  an  example,  in  Fig.  13,  we  have  shown  the  radar  cross  section 
as  a  function  of  K  for  the  Gaussian  distribution  for  0  *  %.  The  cross 
section  was  computed  from  Eqs.  (9)  and  (24)  for  K  s  and  from  the 
Born  approximation  for  K  5  For  K  =  %,  the  two  methods  differed 

by  only  10  per  cent. 


Fig.  13— Radar  cross  section  as  a  function  ofK 
for  fixed  small  positive  j3 
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0  >  1 

For  0  >  1,  we  can  use  Eq.  (24)  to  calculate  cot  11  for  K  <  1 

and  neglect  the  other  phase  shifts.  For  K  >  0,  we  can  use  the  Born 

approximation  to  the  differential  cross  section.  If  0  is  large,  then 

2 

the  geometrical -optics  approximation  can  be  applied  for  1  <  K  <0. 

For  all  other  values  of  K  >  1,  we  can  calculate  the  individual  phase 
shifts  by  the  WKB  approximation  when  the  phase  shifts  are  large,  and 
by  the  small-perturbation  (Born)  approximation  when  the  phase  shifts 
are  small. 

As  an  example,  we  have  calculated  the  radar  cross  section  as  a 

function  of  K  for  the  Gaussian  distribution  for  0  *  18.  We  used 

Eqs.  (9)  and  (24)  for  K  <  1  and  the  WKB  method  for  K  ^  1.  The  results 

are  shown  in  Fig.  14.  In  the  same  diagram,  we  have  also  plotted  the 

geometrical  optics  radar  cross  section  for  \  *  K  s  4.  This  very  useful 

approximation  is  apparently  valid  for  K  as  small  as  unity,  provided 
2 

that  0/K  is  sufficiently  large. 
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Fig.  14— Radar  cross  section  as  a  function  of  K 
for  fixed  large  positive  (3 


